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ABSTRACT 

We present a comprehensive study of the emission spectra from accreting 
sources. We use our new reflection code to compute the reflected spectra from an 
accretion disk illuminated by X-rays. This set of models covers different values of 
ionization parameter, solar iron abundance and photon index for the illuminating 
spectrum. These models also include the most complete and recent atomic data 
for the inner-shell of the iron and oxygen isonuclear sequences. We concentrate 
our analysis to the 2 — 10 keV energy region, and in particular to the iron K-shell 
emission lines. We show the dependency of the equivalent width (EW) of the Fe 
Ka with the ionization parameter. The maximum value of the EW is ~ 800 eV 
for models with log £ ~ 1.5, and decreases monotonically as £ increases. For 
lower values of £ the Fe Ka EW decreases to a minimum near log £ ~ 0.8. We 
produce simulated CCD observations based on our reflection models. For low 
ionized, reflection dominated cases, the 2 — 10 keV energy region shows a very 
broad, curving continuum that cannot be represented by a simple power-law. 
We show that in addition to the Fe K-shell emission, there are other prominent 
features such as the Si and S La lines, a blend of Ar vm-xi lines, and the Ca x 
Ka line. In some cases the S xv blends with the He-like Si RRC producing a 
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broad feature that cannot be reproduced by a simple Gaussian profile. This could 
be used as a signature of reflection. 


1. Introduction 

Accreting systems are observed to emit copious radiation in the X-ray energy range 
which suggests emission from the innermost regions of an accretion disk. Analysis of the 
X-ray spectra is crucial to study the complex mixture of emitting and absorbing components 
in the circumnuclear regions of these systems. However, there are only few observables that 
provide indication of the existence of an accretion disk. These need to be understood in order 
to correctly interpret the physics of these systems. In a general sense, the observed radiation 
from accreting sources can be divided into: a thermal component, in the form of a black body 
emitted at the surface of the disk with typical temperatures of ~ 0. 01 — 10 keV (for the mass 
range ~ 10 8 — 10 M 0 ); a coronal component in the form of a power law covering energies up 
to ~ 100 keV, believed to arise from inverse Compton scattering in a hot gas that lies above 
the disk; and a reflected component, resulting from the interaction of some of the coronal 
X-ray photons and the optically thick material of the disk. In the reflected component, the 
most prominent feature is the iron Kcc emission line at ~ 6.4 keV, produced by transitions 
of electrons between the Is and 2 p atomic orbitals. These are ubiquitous in the spectra of 
accreting sources (Pounds et al. 1990; Nandra &; Pounds 1994; Miller 2007). Other reflection 
signatures are the so called Compton shoulder (next to the Fe K-line), and the Compton 
hump (above ~ 10 keV), produced by the down-scattering of high energy photons by cold 
electrons. 

Much theoretical effort has gone into studies of X-ray illuminated disks over the past few 
decades. Most models assume that the gas density is constant with depth (Done et al. 1992; 
Ross & Fabian 1993; Matt et al. 1993; Czerny & Zycki 1994; Krolik et al. 1994; Magdziarz 
& Zdziarski 1995; Ross et al. 1996; Matt et al. 1996; Poutanen et al. 1996; Blackman 1999). 
Although constant density models may be appropriate for radiation-pressure dominated 
disks, other studies have shown significant differences when the gas density is properly solved 
via hydrostatic equilibrium (Rozanska & Czerny 1996; Nayakshin et al. 2000; Nayakshin &; 
Kallman 2001; Ballantyne et al. 2001; Dumont et al. 2002; Ross & Fabian 2007). Recently we 
have developed a new model for the reflected spectra from illuminated accretion disks called 
xillver (Garcia & Kallman 2010). Although our code is similar in its principal assumptions 
to previous models, xillver includes the most recent and complete atomic data for K-shell 
of all relevant ions (Kallman et al. 2004; Garcia et al. 2005; Palmeri et al. 2008; Garcia et al. 
2009; Palmeri et al. 2010). This has a dramatic impact on the predicted spectra, in particular 
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the K a emission from iron. With this model we can study the effects of incident X-rays on 
the surface of the accretion disk by solving simultaneously the equations of radiative transfer 
and ionization equilibrium over a large range of column densities. Plane-parallel geometry 
and azimuthal symmetry are assumed, such that each calculation corresponds to an annular 
ring at a given distance from the source of X-rays. The redistribution of photons due to 
Compton down-scattering is included by using a Gaussian approximation for the Compton 
kernel. With xillver we are able to solve the reflection problem with great detail, i.e., with 
very high energy, spatial and angular resolution. 

In this paper we present a systematic analysis of our models for reflected spectra from 
X-ray illuminated accretion disks. We show how the most relevant atomic features in the 
spectra depend on the assumed properties of the irradiated gas. We pay particular attention 
to the Fe K-shell emission lines, and we quantify its strength in terms of the equivalent widths 
predicted by our models. These models are also used to produce faked CCD spectra in order 
to simulate the effects of the instrumental response and limited spectral resolution. These 
results will be helpful diagnostic tools in the interpretation of accreting sources observations. 

In the next Section we describe briefly the numerical methods used in our reflection 
code. In Section 3 we present the results of our analysis on the simulated spectra, as 
well as comparisons with observations from Seyfert galaxies and X-ray binaries. The main 
conclusions are presented in Section 4 


2. The Reflection Model 


In order to calculate the reflected spectra from X-ray illuminated accretion disks we 
make use of our reflection code xillver. The details of the calculations are fully described 
in Garcia &; Kallman (2010), thus here we just review the main aspects. 

The description of the interaction of the radiation with the gas in the illuminated slab 
requires the solution of the transfer equation: 


2 2 fw ^ 
ji u ( E,t ) 


dr 2 


+ /i 2 uj(E, t ) 




dr 
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where u(/ji,E,t) is the average intensity of the radiation field for a given cosine of the 
angle with respect to the normal //, energy E. and position in the slab, specified by the 
Thomson optical depth dr = —a^dz — — aTn e dz , where ot is the Thomson cross section 
(= 6.65 x 10 -25 cm 2 ), and n e is the electron number density. The quantity ui is defined as 
the ratio of the Thomson scattering coefficient a T to the total opacity y(r, E). The second 
term in the right-hand side of Equation 1 is the source function, which is defined as the ratio 
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of the total emissivity over the total opacity of the gas, taking into account both scattering 
and absorption processes: 
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where ctk n (E) is the Klein-Nishina scattering coefficient, j(E,r ) is the thermal continuum 
plus lines emissivity, and J C (E, r) is the Comptonized mean intensity, given by the Gaussian 
convolution 
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The Gaussian is centered at E c — E'{ 1 + 49 — E'/m e c 2 ), where m e is the electron mass, c is 
the speed of light, and 9 — kT/m e c 2 is the dimensionless temperature. The energy dispersion 
is given by a = E' [29 + \{E' jm e (?f ] 1/2 . 


The solution of the system is completed by imposing two boundary conditions. At the 
top of the slab (r = 0), we specify the incoming radiation field incident at a given angle [io 

by 
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where F X (E) is the net flux of the illuminating radiation incident at the top of the slab. 
At the inner boundary (r = T max ), we specify the outgoing radiation field to be equal to a 
blackbody with the expected temperature for the disk: 


u(T max ,H,E)n 
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where T disk can be defined using the Shakura & Sunyaev (1973) formulae. Since these models 
are calculated under the assumption of constant density, we use the common definition of 
the ionization parameter (Tarter et al. 1969) to characterize each case: 


^ _ 47 rF x 
n e 

where n e is the electron gas density, and F x is the net illuminating flux integrated in the 
1 — 1000 Ry energy range. The solution of the system is found by forward elimination and 
back substitution. A full transfer solution must be achieved iteratively in order to self- 
consistently treat the scattering process. This procedure requires ~ T^ ax iterations (lambda 
iterations) for convergence. 



Given the solution for the radiation field at each point in the atmosphere, we use the 
photoionization code XSTAR (Kallman k, Bautista 2001) to determine the state of the gas 
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at each point of the gas. The state of the gas is defined by its temperature and the level 
populations of the ions. The relative abundances of the ions of a given element and the level 
populations are found by solving the ionization equilibrium equations under the assumption 
of local balance, subject to the constraint of particle number conservation for each element. 
XSTAR calculates level populations, temperature, the opacity y(i?, r) and the total emissivity 
j(E, t) of the gas assuming that all the physical processes are in steady state and imposing 
radiative equilibrium. 

The XSTAR atomic database collects recent data from many sources including CHIANTI 
(Landi & Phillips 2006), ADAS (Summers 2004), NIST (Ralchenko et al. 2008), TOPbase 
(Cunto et al. 1993) and the IRON project (Hummer et al. 1993). The database is described 
in detail by Bautista & Kallman (2001). Additionally, the atomic data associated with the K- 
shell of the Fe ions incorporated in the current version of XSTAR has been recently calculated 
and represents the most accurate and complete set available to the present. A compilation 
of these results and a careful study of their impact on the photoionization models can be 
found in Kallman et al. (2004). Moreover, XSTAR also includes the atomic data relevant to 
the photoabsorption near the K edge of all oxygen (Garcia et al. 2005), and nitrogen ions 
(Garcia et al. 2009). 

Finally, all the calculations presented in this paper were carried out over a large optical 
depth ( T max — 10), considering high resolution spectra with an energy grid of at least 5 x 10 3 
points (1Z — E/AE ~ 350), 200 spatial zones, and 50 angles to account for anisotropy of 
the radiation field. The simulations do not take into account the dynamics of the system. 
Input parameters common to all these models are: the electron gas density n e — 10 15 cm -3 , 
photon index of the incident radiation T — 2, mass of the central object M — 10 8 M 0 , 
distance from the central object R — 7 R s , and the mass accretion rate M — 1.6 x 10 ~ 3 M E dd, 
where R s — 2 GM/c 2 is the Schwarzschild radius and M E dd is the accretion rate at the 
Eddington limit. This particular set of parameters yields a disk effective temperature of 
Tdisk — 2-8 x 10 4 °K. The models presented here cover 10 different values for the illumination 
flux F x = 5 x 10 14 — 5 x 10 17 erg cm -2 s -1 , which corresponds to ionization parameters of 
log £ = 0.8, 1.1, 1.5, 1.8, 2.1, ..., 3.8. It is important to notice that these simulations are not 
constrained to a particular geometry for the illumination. This is because the flux F x of the 
illuminating radiation is defined at the surface of the slab, regardless of the position of the 
X-ray source. Therefore, a set of models with different values of the ionization parameter can 
be used to construct a particular geometry, by defining the luminosity and location of the 
source. We do not consider models with ionization parameters lower than ~ 4ir (log £ = 0.8), 
since our code is optimized for calculations of medium to high ionization. The upper limit is 
set to log £ = 3.8, because at such a high illumination the gas is almost completely ionized 
over a large depth. 
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3. Results 


3.1. Equivalent Widths 


The equivalent width of a line provides a quantitative measure of the strength of the 
spectral profile, both in emission or absorption. It is defined by the well known formula 


EW = 


pE 

Je 1o 


E high ( F (£) _ F C (E)) 


F C (E) 


dE, 


(7) 


where F(E) is the total flux and F C (E) is the total flux of the continuum under the line. 
The integration is performed over the energy range where the spectral feature takes place, 
between Ei ow and E high . Usually, there are uncertainties in the determination of the intrinsic 
continuum which affects the knowledge of the integration region and the actual value of the 
equivalent width itself. Nonetheless, one can approximate its calculation by defining an 
energy region in which one knows only the spectral feature of interest appears, and where a 
local continuum can be defined. 


In this paper we make use of the equivalent width as a measure of the strength of several 
features in the X-ray spectra reflected from illuminated accretion disk. As an example, in 
Figure 1 we show the reflected spectra from 3 different models calculated for log £ = 0.8, 1.8 
and 2.8 in the 4-9 keV energy region, where the only atomic features are due to inner-shell 
transitions from Fe ions. The spectra (flux vs. energy) are shown as solid lines. Vertical 
dotted lines are placed at 5.5 keV and 7 keV, defining a particular integration region. The 
continuum is defined as a straight line that passes through these two points in the spectra, 
which is shown as dashed lines. It is clear from the Figure that the resulting continuum does 
a good job reproducing the local continuum and that it is not superimposed over any part 
of the emission profile. The upper limit Ehigh = 7 keV was chosen such that only emission 
from Ka transitions of Fe are taken into account, neglecting any K/3 emission (which occurs 
at energies above 7 keV). The lower limit, however, is more arbitrary. It needs to be chosen 
such that all the line profile is included in the integration. If Ei ow is set too large, part 
of the line emission can fall outside the range, especially for cases with large values of 
the ionization parameter, where Comptonization smears the line profile due to the down- 
scattering. However, if E iow is set to a small value, all the line emission is taken into account 
but the local continuum is modified to a point that could either under- or over-estimate the 
strength of the line. Therefore, we have repeated the calculation of the equivalent widths 
for Ei ow — 5,5.5 and 6 keV, in order to evaluate the sensitivity of the results on this lower 
limit. 


Figure 2 shows the resulting equivalent widths as a function of the ionization parameter, 
calculated for the Fe Kq; emission profile in 10 of our models. Circles connected with solid 
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lines correspond to the integration range 5-7 keV, squares connected with dashed lines to 
5.5-7 keV, and triangles connected with short-dashed lines to the integration performed 
over the 6-7 keV energy range. It is clear from the Figure that the equivalent widths are 
almost unaffected by the variations of the lower boundary of the integration range, with 
the exception of those models with 1.5 <log £ < 2.5. These are the most sensitive cases 
probably because of the complexity of the iron K emission. For lower values of the ionization 
parameter, the emission mainly occurs at 6.4 keV due to mostly neutral Fe ions. Higher 
values of £ means that the gas is highly ionized and thus mostly He- and H-like Fe ions are 
responsible for the emission at ~ 6.9 keV. In both cases the line profile is simple in the sense 
that the emission is concentrated at one particular energy. In between, emission from many 
different ions takes place at the same time, creating a more complex spectral feature, as can 
be seen in the spectrum for log £ = 1.8 in Figure 1. Since there is no clear reason to choose 
one of these values of Ei ow , we choose the intermediate one (E iow = 5.5 keV). We consider 
the uncertainties in the equivalent widths to be of order of the differences between the values 
shown in Figure 2 (< 100 eV). 

Figure 3 shows a comparison of the Fe Ko: equivalent widths predicted by our models 
(circles connected with dashed lines), and those predicted by the models included in REFLION 
(triangles connected with solid lines, Ross &: Fabian 2005). The calculation of the equivalent 
widths is the same in both models, with the integration performed in the 5.5-7 keV energy 
range. There are important differences to notice in this comparison. It is convenient to make 
the distinction between two regions in the plot; the region for models with log £ > 1.5 and 
the one for models with log £ < 1.5. In the high ionization region, both sets of models show 
a similar behavior, the equivalent widths decrease monotonically as the ionization parameter 
increases, resembling the Baldwin effect for X-rays. Nevertheless, in this region all our models 
systematically predict stronger Fe Ka: emission. For those models in the low ionization region 
(log £ < 1.5), the differences are significantly larger. In fact, reflion models predict that 
the equivalent widths keep growing as the ionization parameter decreases, while our models 
show the opposite trend. This turn over in the values of the equivalent widths can be 
understood by looking the ionization balance in detail. Kallman et al. (2004) performed 
similar calculations in photoionized models using the same atomic data used in our models. 
Using a set of XSTAR models of thin spherical shells illuminated by the same power law 
spectrum they calculated the ratio of the emissivity per particle for K line production as 
a function of the ionization parameter. Their Figure 7 shows the contribution of each ion 
of iron summing over the K lines for each one. The overall behavior of these models is 
quite similar to the ones presented in this paper. The combined line emissivities have a 
minimum just below log £ ~ 1 (roughly where our calculations start), and is mainly due to 
a combination of Fe xiv-xvi. As the ionization increases the also the emissivities increase, 
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with maximum values that peak around log £ ~ 2, where the emission is dominated by 
Fe xvill and Fe xix. For log £ > 2.5 the overall line emission decreases rapidly as the iron 
ions become more ionized. This clearly resembles the general trend exhibit by the equivalent 
widths derived from our models, as is shown in Figure 3. It is worthwhile to mention that the 
equivalent widths will increase again for log £ < 1, as can be seen in Figure 7 of Kallman et al. 
(2004), but those values of the ionization parameter are out of the range of the calculations 
presented here. The sensitivity of the iron Ka equivalent width to the ionization parameter 
represents an important potential diagnostic of photoionized plasmas, given its ubiquity in 
the observed spectra from accreting sources. 

Another important diagnostic can be achieved through the dependence of the equivalent 
width of the Ka line of iron on its abundance with respect to the solar value. This can be seen 
in Figure 4, where the equivalent width resulting from models with different iron abundances 
is plotted. In the Figure, connecting lines correspond to a particular ionization parameter 
(shown next to each curve). For each value of £, three reflected spectra were calculated 
assuming 0.2, 1 and 10 times the solar abundance of iron. The resulting equivalent widths 
are shown as filled circles for each case. It is clear from Figure 4 that the dependence of 
the Fe Ka equivalent width on the ionization parameter is almost the same for different 
abundances of iron, i.e., the shape of the curve in Figure 3 will remain unchanged if we use 
models with iron over or under abundant, and only the actual values of the equivalent widths 
will show a variation. This is a consequence of the fact that the Fe Ka equivalent width 
grows linearly with the iron abundance, which is also evident from the Figure. To illustrate 
this, we include the function y(x ) = lOx with a dashed line in Figure 4, which shows a slope 
similar to all the other curves. 


3.2. Simulated observations 

The analysis presented previously provide important tools for the interpretation of the 
X-ray spectra from accreting sources by the use of simple quantities that can be easily de- 
rived from the observed spectrum. Nevertheless, the models computed here do not necessarily 
represent real observed data. First, the spectral resolution used in the models ( R ~ 350), 
although comparable to grating spectra, is far superior to what is achieved in CCD ob- 
servations. Second, the models do not suffer from loss of information due to the natural 
restrictions involved in real observations, such as the instrument effective area which affects 
the total number of photons collected. Finally, a real observation is, in general, a mixture 
of several components and not exclusively represented by a reflected spectrum. In accreting 
sources such as active galactic nuclei (AGN) and galactic black holes (GBHs), an important 
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component is the contribution of the source of X-rays that illuminates the accretion disk. In 
fact, it is not clear which fraction of the direct over the reflected component can be present 
in a given observation, and this proportion is most likely to depend on the geometry and 
orientation of the system. 

In order to account for some of these effects, we have produced simulated CCD observa- 
tions using our models as the source of the reflected spectra. We have used the fake it task 
in XSPEC in combination with the FI XIS response matrices from Suzaku to simulate 100 
ksec observations, assuming a source flux of 3 x 10 -10 erg cm -2 s -1 in the 2-10 keV energy 
band. This flux corresponds to 10 mCrab, which is a typical value for bright AGN. We then 
try to fit the resulting simulated observation with a simple phenomenological model, i.e., 
the basic continuum is fitted with a power law, and all the relevant features with Gaussians. 
An example of such a fit is shown in Figure 5. In the upper panel, data points with error 
bars represent the simulated data using a reflection model with log £ = 0.8, while the solid 
line shows the best-fit. The lower panel shows the theoretical model used with all its com- 
ponents. The continuum is fitted with a power law (thick dashed line), and several emission 
profiles with Gaussians (thin dotted lines) centered at 2.39, 3.06, 3.8, 6.38 and 7.03 keV, cor- 
responding to emission from S Vl-x, Ar VIII-XI, Ca X, Fe x-xil Ka and Fe K/3, respectively. 
Finally, a very broad (a — 1.54 keV) Gaussian profile centered at 6.1 keV was included in 
order to properly fit the continuum at lower energies. Clearly, a simple power law cannot 
describe the continuum for models with low ionization, given the significant modification of 
the original power law continuum due to the large values of the photoelectric opacity. This 
very broad curving continuum is required to fit the reflected spectra with log £ ~ 1.1 or 
lower. Therefore, this particular profile could be used as a signature of strong reflection from 
low ionized atmospheres. 


We have then produced simulated observations using all the models considered previ- 
ously, i.e., using the 10 reflected spectra calculated for log £ = 0.8, 1.1, 1.5, 1.8, ..., 3.8. By 
including the Suzaku response matrices we can simulate the effects of the effective area from 
the instrument, its sensitivity and energy resolution. However, many X-ray observations 
from accreting sources reveal that the direct, unprocessed radiation as a power-law is de- 
tected simultaneously along with the reflected component. Therefore, we have produced 
additional simulated observations with a source spectrum resulting from the combination of 
both the reflected and the illuminating components. In order to quantify the proportion of 
one component to the other we use a (3 parameter, defined as the ratio of the incident over 
the reflected flux, 


P = 



( 8 ) 


where both fluxes are defined in the 2-10 keV energy range. Note that this quantity differs 
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from the reflection fraction R — Q/27T, where Q is the solid angle sustained by the reflector. 
The relation between the reflection fraction and our (3 parameter is roughly R rs 1//3; such 
that a pure reflection case corresponds to (3 — 0 and R ~^> 1. while in the case of no 
reflection (3 1 and R — 0. Because our models do not carry any information about the 

geometry of the reflector, it is not appropriate to derive a reflection fraction R for them. 
Instead, we use the (3 parameter to quantify the dilution of the spectral profiles by the direct 
power-law spectrum, in order to compare derived spectral quantities with the those from 
real observations. 

Figures 6, 7 and 8 show the simulated CCD observations using our models with log £ = 
0.8, 1.8 and 2.8, respectively. In each figure, the data points in the upper panel show the 
simulated data for three different values of the (3 parameter, namely, from bottom to top, 
(3 — 0, 0.5 and 2. The solid lines are the best-fit for each case, using the same kind of simple 
models described before. The dots in the lower panel represent the residuals in units of a 
for all the fits. These figures show the effect of the instrument on the observed spectra. In 
particular, by comparing the curve for (3 — 0 in Figure 6 with the corresponding model in 
Figure 1, it is clear how the narrow Fe K a and K(3 lines in the log £ = 0.8 model become 
much broader features in the simulated observation. As mentioned before, emission from S, 
Ar and Ca ions are the most prominent features at energies < 6 keV. It is also important 
to notice from this Figure the dilution of the spectral profiles due to the inclusion of the 
direct power- law spectra (i.e., the curves for [3 — 0.5 and (3 — 2). The iron Kck at 6.4 keV 
and the sulfur blend at ~ 2.39 keV emission profiles are clearly detectable, even for (3 = 2 
(i.e., twice as much direct flux as the reflected flux in the 2-10 keV band). However, the 
weaker emission lines become almost undetectable for (3 > 2. Most important, not only 
the atomic features are modified but the continuum as well. The broad Gaussian profile 
required to fit the continuum in combination with the power-law is no longer required for 
f3 — 0.5 or greater. This is important, since this means that this feature is characteristic of 
low ionization, but also appears exclusively for the pure reflection cases ((3 — 0). 

Figures 7 and 8 show essentially the same situation. In Figure 7, which corresponds to 
the simulation generated with log£ = 1.8, it is important to notice the presence of the S xv 
Kq; emission line at ~ 2.45 keV, and the S xvi Ka plus the radiative recombination continua 
(RRC) at ~ 2.59 keV, which are detectable even in the (3 — 2 simulated spectrum. The iron 
K (3 emission line becomes undetectable for (3 > 1. Figure 8 shows the log£ = 2.8 simulated 
spectra. In this case, the iron emission is the only clear and strong feature. However, the iron 
K-shell emission can no longer be well represented by a single Gaussian profile. Instead, a 
broad (a — 1.26 keV) Gaussian centered at ~ 6.5 keV accounts for the Comptonized profile, 
while a narrow one centered at ~ 6.6 keV accounts for the discrete emission lines that can 
be seen in the original spectrum of Figure 1. The iron emission can be well fitted even in the 
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(3 — 2 case, but not for cases with a larger contribution of the direct power-law spectrum. 

In Table 1 we show a compilation of the strongest features detectable in the simulated 
CCD spectra for the pure reflection case (/ 3 — 0). We show the line centroid energy and 
the equivalent with for each feature. Three of the features are blends of several components, 
such as the Si Vl-x around 2.39 keV, the blend of the S XI Ko; line and the Si xiv RRC 
around 2.59 keV, i and the Ar viii-xi blend around 3.3 keV. 


3.3. Comparison with observations 

In order to test the validity of the results obtained with the simulated CCD spectra, 
we have compared our estimates with those derived from real observations. To do this we 
have concentrated our attention on the equivalent widths and centroid line energies of the 
iron K a. emission profile, since this is the strongest emission present in the spectra from 
most accreting sources. Because the resulting equivalent width is very sensitive to small 
changes in the fitted continuum, we decided to apply a systematic procedure to calculate the 
line equivalent width, similar to the one described in Section 3.1. First, we have generated 
simulated CCD spectra for the 10 different values of the ionization parameter considered 
in this paper (log £ = 0.8, 1.1, 1.5, 1.8, ..., 3.8), and for 5 different values of the (3 parameter 
from Equation 8 ((3 — 0, 1, 2, 5, 10), which covers cases from pure reflection cases up to cases 
dominated by the direct power-law (no reflection). Then, we fit a simple model to each 
spectrum, consisting of a power-law for the local continuum and one Gaussian profile for the 
Fe K-shell emission. The fit is performed in a small energy range concentrated around the 
iron lines, ignoring any other channels. We choose the range to be 5. 7-7.3 keV for spectra 
with 0.8 > log £ > 2.8, and 3-10 keV for those with larger ionization parameters, since the 
Fe K-line is much broader due the Comptonization. 

In Figure 9, the line centroid energies and equivalent widths resulting from all our 
simulated spectra are compared with the same quantities derived from spectra of Seyfert 
(Sy) galaxies and low mass X-ray binaries (LMXB). In the Figure, filled circles correspond 
to Sy 1 and 1.2 galaxies, asterisks to Sy 1.5, and filled triangles to Sy 1.8, 1.9 and 2; all 
taken from the Winter et al. (2009) compilation. Filled squares correspond to the LMXB 
analyzed by Ng et al. (2010). The error bars are included as dashed lines. Dots connected 
with solid lines are the values predicted by our reflection models. Each line corresponds to 
one particular ratio of the direct over the reflected component, or (3 parameter. From top to 
bottom, each curve corresponds to (3 — 0, 1, 2, 5 and 10. Error bars are also included as solid 
lines. Along each line, the points from left to right correspond to the models with increasing 
values of the ionization parameter £. 
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There are important pieces of information contained in Figure 9. The observed data 
shows the clear distinction between Sy galaxies and LMXB with respect to their iron emis- 
sion. Sy galaxies seem to only show the cold iron line at ~ 6.4 keV, while in the LMXB 
spectra the hot line around 6.7 keV is detected. This yields a marked difference in the ion- 
ization of the emitting gas between these two types of sources. According to our models, 
the cold region of the plot where Sy galaxies appear (for line energies lower than 6.5 keV), 
correspond to values of log £ < 2, while the hot region (line energies larger than 6.5 keV), 
where LMXB show up, covers larger values of ionization, up to log £ = 4. In the cold region, 
Sy galaxies equivalent widths cover a large range, from 10 to 1000 eV. Despite the large error 
bars in many of the values, equivalent widths larger than 200 eV seem to correspond almost 
exclusively to Sy 1.5-2 galaxies. Our models predict that the largest values correspond to re- 
flection dominated cases, with (3 — 0,1 and 2. Lower equivalent widths (less than 200 keV), 
seem to be common to all types of Sy galaxies. However, our models require significant 
contribution from the direct power-law illuminating spectra to dilute the line profiles and 
reproduce the lower values of the equivalent widths. Accordingly, we need to produce models 
with (3 — 10 to explain the equivalent widths lower than 100 eV (i.e., 10 times more direct 
power- law flux than the reflected one, or R — > 0). Also, it is important to notice that for 
the lower ionization models the resulting line energies are much closer to each other than for 
the high ionization ones. At the same time, the data for Seyfert galaxies in Figure 9 shows 
a large concentration of points at 6.4 keV. This is probably due to the fact that sometimes 
observers fix the line centroid energy to a known value given the poor statistics of the data. 
This imposes difficulties in the estimation of the ionization parameter when comparing the 
data values to the ones predicted by the models. Therefore, accurate line energies are needed 
from the observations in order to improve the diagnostics. 

The LMXB observations are more scattered in the hot region of the Figure 9 (line 
energies larger than 6.5 keV). Both the models the the observed data predict equivalent 
widths as small as ~ 10 eV, but no larger than ~ 200 eV. The error bars in the values 
resulting from the models are much larger than in the low ionization cases, mainly because 
the line feature is weaker as the ionization increases, since the emission is mostly due to 
H- and He-like iron. The line profile is even weaker as the (3 parameter increases, due to 
the direct power-law dilution. Therefore, in the cases lor (3 — 1 and 2 (third and fourth 
curves from top to bottom), the line cannot be detected in the spectra simulated with 
highest ionization parameter (log £ = 3.8). The same occurs for the simulated spectra with 
log £ = 3.5 in the bottom curve (/ 3 — 10). The ranges of the ionization and the (3 parameters 
covered in this paper are sufficient to well represent the observed data for both Seyfert 
galaxies and LMXB. 

The fact that the largest equivalent widths predicted by our reflection dominated models 
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correspond mostly to the values found in Sy 1.5-2 galaxies agrees with the general picture 
that in these sources the direct component of X-ray radiation is being obscured, and only the 
reflected or reprocessed radiation is observed. Following this idea, we repeated the compar- 
ison of our models with a larger and more recent compilation of AGN Seyfert galaxies done 
by Fukazawa et al. (2010). These authors analyzed high-quality Suzaku data for 88 Seyfert 
galaxies, and derived line energies, equivalent widths, and hydrogen column densities among 
other quantities. In Figure 10 we show the equivalent width against the Fe Ka line energy 
derived from these observations. For comparison we use the same scale used in Figure 9, 
but in this case we use different symbols to show different values of the hydrogen column 
density JV H . In the plot, filled circles represent galaxies with log 7V H <21 cm -2 , asterisks 
are galaxies with 21 < jV H < 23 cm -2 , and filled triangles are those with log 7V H > 23 cm -2 . 
As in Figure 9, filled squares are the values for the LMXB from Ng et al. (2010), and dots 
connected with solid lines are the values predicted by our models. This particular sample 
has the advantage of being more accurate (lower error bars), and is more consistent since 
all the data was obtained with the same instrument and analysis technique. However, the 
information derived from this comparison is consistent with the Winter et al. (2009) compila- 
tion. Equivalent widths larger than 200 eV are only observed in sources with large hydrogen 
column densities (log N n > 23 cm -2 ), which corresponds to the reflection-dominated models 
(0</3< 2). Lower equivalent widths are observed in sources with all values of Nu- The 
majority of sources with low hydrogen column densities (filled circles) have equivalent widths 
lower than 100 eV, consistent with models with large contribution of the direct power-law 
(ft ~ 10). The agreement between these observations and our models is in general very good. 
However, there are at least five galaxies with equivalent widths larger than the maximum 
values predicted by our models (~ 1 keV) with solar abundances. However, larger equivalent 
widths can be easily achieved by increasing the iron abundance, as shown in Figure 4. 


4. Conclusions 

In this paper we have presented a comprehensive study of the emission spectra from 
accreting sources. Using our new reflection modeling code (Garcia & Kallman 2010), we 
have computed the reflected spectra from an X-ray illuminated accretion disk. We have 
concentrated our analysis to the iron K-shell emission lines, although other spectral features 
are also discussed. Our models predict equivalent widths for the Fe Ka emission of ~ 10 eV 
for high ionization parameters (log £ ~ 4), and maximum values of ~ 800 eV for models 
with log £ ~ 1.5. For lower values of the ionization parameter the equivalent widths decrease 
to a minimum near log £ ~ 1, contrary to what other models predict. These differences are 
due to the atomic data used by each simulation. We have shown that the behavior of the Fe 
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Ka equivalent widths with respect to the ionization parameter of the gas is consistent with 
the line emissivities shown by Kallman et al. (2004), where the same atomic data was used. 
Additionally, these equivalent widths display a linear dependency with the iron abundance 
normalized to its solar values. This seems to be true for all the models within the range of 
ionization parameter values considered in this paper. 

Simple analysis of simulated CCD spectra reveals that for low ionized, reflection domi- 
nated cases the 2-10 keV continuum cannot be represented by a simple power-law. Instead, 
a broad Gaussian profile is also required. This type of continuum can be used as a strong 
reflection signature. These simulations also indicate that in addition to the iron K-shell lines, 
the S xv Ka emission line at 2.45 keV is one of the most prominent features in the spectra. 
For cases with log £ ~ 2, this line blends with the He- like silicon RRC providing a broad 
feature that cannot be represented by simple Gaussian profiles. This could also be used as 
another reflection signature while analyzing observations. 

The iron Ka equivalent widths and line centroid energies are in good agreement with 
values reported in the literature for both AGN Seyfert galaxies, and galactic sources such 
as LMXB. In particular, both model and observations show that large equivalent widths 
(> 200 eV) can be achieved in many situations. According to our models, the largest values 
of the equivalent widths correspond to the lower (3 parameters, i.e., reflection dominated 
cases. Only observational data from AGN Sy 1.8-2 galaxies, or those with log iV H > 23 cm -2 
coincide with these large equivalent widths, but only for low values of the ionization param- 
eter (log £ < 2). The majority of the Sy 1-1.2, or those sources with log /V H < 21 cm -2 show 
equivalent widths below 200 eV. In order to reproduce those values with our models, we 
need to increase the contribution of the direct power-law up to ten times the contribution of 
the reflected component. These results agree with the general idea that type 1 and 2 Seyfert 
galaxies only differ on their orientation with respect to the observer. In the case of LMXB, 
both the model and the observed data coincide in equivalent widths of ~ 10 — 200 eV. Since 
all these sources show hotter Fe K-lines (with energies around 6.7— 6.9 keV), they correspond 
to models with values of the ionization parameter larger than those for AGN (2 < log £ < 4). 


This work was supported by a grant from the NASA astrophysics theory program 05- 
ATP05-18. This research has made use of NASA’s Astrophysics Data System. 
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Fig. 1. — Solid curves: emergent spectra from the surface of an accretion disk with constant 
density (n — 10 15 cm -3 ) calculated with our reflection code xillver for three different illu- 
minations. The log of the ionization parameter is shown next to each curve. The curves are 
rescaled for clarity. Dashed lines: local continuum used in the calculation of the equivalent 
widths. Dotted vertical lines are placed at 5.5 and 7 keV. 
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Fig. 2. — Equivalent width of the Fe Ka line as a function of the log of the ionization 
parameter, resulting from three different energy ranges of integration. Connected circles, 
squares and triangles correspond to Et ow — 5, 5.5 and 6 keV, respectively. E^igh is set equal 
to 7 keV for all the cases. 
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Fig. 3. — Equivalent width of the Fe Ka line as a function of the log of the ionization 
parameter. Connecting circles are the values predicted by our reflection code xillver, 
while triangles correspond to the values predicted by reflion (Ross &; Fabian 2005). The 
integration is performed over the 5.5-7 keV energy range in both models. 
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Fig. 4. — Fe Kq; equivalent width as a function of the iron abundance (with respect to solar). 
Connecting circles are the values predicted by xillver for different illuminations. The value 
of log £ is shown next to each curve. The dashed line represents the function y(x) — lOx, to 
exemplify the linear dependency of these quantities. 
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Energy (keV) 

Fig. 5. — Upper panel: data points are the simulated CCD observation produced with our 
model for log £ = 0.8. The solid line is the best-fit using a simple phenomenological model. 
Error bars are shown for all the data points. Lower panel: the theoretical model used in the 
fit is shown with the solid line, while all the individual components are shown with dashed 
curves for the Gaussian profiles, and with a thick-dashed line for the power-law. 
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Fig. 6. — Upper panel: data points show the simulated CCD spectra produced with our 
model for log £ = 0.8, for (3 — 0, 0.5 and 2, from bottom to top. The solid lines are the 
best-fits for each case. Lower panel: residuals in units of a for all the fits. Error bars are 
shown for all the data points. 
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Fig. 7. — Upper panel: data points show the simulated CCD spectra produced with our 
model for log £ = 1.8, for f3 — 0, 0.5 and 2, from bottom to top. The solid lines are the 
best-fits for each case. Lower panel: residuals in units of a for all the fits. Error bars are 
shown for all the data points. 
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Fig. 8. — Upper panel: data points show the simulated CCD spectra produced with our 
model for log £ = 2.8, for (3 — 0, 0.5 and 2, from bottom to top. The solid lines are the 
best-fits for each case. Lower panel: residuals in units of a for all the fits. Error bars are 
shown for all the data points. 
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Fig. 9. — Equivalent width versus centroid energy of the Fe Ka line. Filled circles, asterisks 
and filled triangles correspond to AGN Seyfert 1-1.2, 1.5 and 1.8-2 galaxies, respectively 
(Winter et al. 2009). Filled squares correspond to LMXB (Ng et al. 2010). Dots connected 
with solid lines are the values predicted by our models with 0.8 < log £ < 3.8 (from left to 
right). From top to bottom, each curve corresponds to /3 = 0, 1, 2, 5 and 10. Error bars are 
shown for both data and model points. 
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Fig. 10. — Equivalent width versus centroid energy of the Fe Ka line. Filled circles, asterisks 
and filled triangles correspond to AGN Seyfert galaxies with log Nu < 21, 21 < log Nr < 23, 
and log Nu > 23 cm -2 , respectively (Fukazawa et al. 2010). Filled squares correspond 
to LMXB (Ng et al. 2010). Dots connected with solid lines are the values predicted by 
our models with 0.8 < log £ < 3.8 (from left to right). From top to bottom, each curve 
corresponds to /? = 0, 1, 2, 5 and 10. Error bars are shown for both data and model points. 
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Table 1. Strongest features in the simulated CCD spectra 


log £ 

(keV) 

0.8 

B 

B 

1.8 

2.1 

2.5 

2.8 

3.1 

3.5 

3.8 

Si xiv K a 

Eline 




1.98 

1.93 

2.01 

2.01 





EW 




0.19 

0.26 

0.04 

0.02 




Si Vl-X 

^line 

2.39 

2.39 










EW 

6.12 

1.03 









S xv Ka 

Eline 



2.46 

2.45 

2.45 

2.44 





+ Si xm RRC 

EW 



0.36 

0.29 

0.17 

0.03 





S xi Ka 

Inline 



2.58 

2.59 

2.59 

2.59 

2.57 




+ Si xiv RRC 

EW 



0.05 

0.09 

0.10 

0.04 

0.03 




Ar VIII- XI 

Inline 

3.06 

3.07 

3.29 

3.29 

3.34 







EW 

0.63 

0.10 

0.07 

0.11 

0.23 






Ca x 

Eline 

3.81 

3.83 

3.84 









EW 

0.43 

0.14 

0.03 








Fe Ka 

^ line 

6.38 

6.38 

6.39 

6.41 

6.46 

6.61 

6.61 

6.65 

6.74 

6.92 


EW 

0.38 

0.84 

1.07 

0.87 

0.65 

0.42 

0.37 

0.39 

0.16 

0.04 

Fe K/3 

line 

7.03 

7.01 

6.96 

6.97 

6.97 













